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ABSTRACT 

In the ACDM cosmology dark matter halos grow primarily through the accretion of 
smaller halos. Much of the mass in a halo of 10 14 Mq comes in through accretion 
of ~ 10 13 Mq halos. If each such halo hosted one luminous red galaxy (LRG) then 
the accretion of so many halos is at odds with the observed number of LRGs in 
clusters unless these accreted LRGs merge or disrupt on relatively short timescales 
(~ 2 Gyr). These timescales are consistent with classical dynamical friction arguments, 
and imply that 2 — 3 LRGs have merged or disrupted within each halo more massive 
than 10 14 Mq by z = 0. The total amount of stellar mass brought into these massive 
halos by z = is consistent with observations once the intracluster light (ICL) is 
included. If disrupted LRGs build up the ICL, then the hierarchical growth of massive 
halos implies that a substantial amount of ICL should also surround satellite LRGs, 
as suggested by recent observations of the Virgo cluster. Finally, we point out that 
these results are entirely consistent with a non-evolving clustering strength and halo 
occupation distribution, and note that observations of the latter in fact support the 
hypothesis that merging/disruption of massive galaxies does indeed take place at late 
times. 



1 INTRODUCTION 



The formation and evolution of massive red galaxies provide 
a critical testing ground for modern theories of galaxy for- 
mation based on hierarchical merging of dark matter halos. 
Ongoing growth of massive halos via mergers is a generic fea- 
ture of hierarchical models, such as cold dark matter (CDM). 
However evidence for the ongoing assembly of massive galax- 
ies is at best inconclusive. Evolution in the galaxy stellar 
mass and luminosity functions at the ma ssive /luminous end 

1 (e . g. iDrorv et al. 



for d e cades (e.g. Cha ndrasckhar 1943; Ostrikcr fc Tremaind 
1 19751 ; lOstriker fe Hausmanl 1 19771 ; iMerritd Il984l h accurate 
merger times of satellite galaxies have been historically hard 
to calculate, and are poorly constrained observationally. Un- 
fortunately, the problem cannot at present by circumvented 
by brute-force simulations due both to severe resolution re- 
quirements and the uncertain effects of b aryon condensation 
on the survival of satellite halos (see e.g.lMoore et al.lll99 9 ; 
Klypin et all 19991 ; IPiemand et al.ll2004l ; lGao et al.ll2004al lb ; 



appears quite modest since z 
Bundv et al.| 120051 ; iBorch et all |2006|: 
Faber et al.l I2OO6I; IWillmer et alj 120061; iBrown et al 



iFontana et al 



2004 



2006 



2007 



Reed et al.ll2005h . In addition, while DF is usually consid 



Caputi et al.l 120061 ; IWake et all 120061 ) though estimates of 



the mer ger rate of massive galaxies present a less consistent 



picture Ivan Dokkuml 120051 : iBell et al.ll2006l ; iMasiedi et"al 



2006; Wh ite et alj|2007i r If, as theory predicts, massive ha- 



los are constantly accreting halos that are themselves hosts 
of massive galaxies, what is the fate of these accreted galax- 
ies? 

Two physical effects can cause satellite galaxies to 'dis- 
appear' from an observational sample. Tidal forces acting 
on a satellite as it orbits in the host halo potential can 
cause it to disrupt. At the same time, dynamical friction 
(DF) causes a satellite to lose energy to the background 
dark matter halo and eventually causes the satellite to sink 
toward the center and merge with the central galaxy of 
the host halo. While such notions, and their relevance to 
the evolution of galaxies within clusters, have been known 



ered in a collisionless medium (such as dark matter), DF 
acting in a collisional medium (such as intracluster gas) is 
stronger (weaker) than in the collisionless c ase for satellite s 
traveling at supersonic (subsonic) speeds l|Ostrikerl Il999l ). 
Observational constraints on the merging timescale of satel- 
lites would hence provide valuable insight into this complex 
dynamical process. 

This paper explores observational constraints on the av- 
erage merging timescale of luminous red galaxies (LRGs). 
We assign LRGs to dark matter halos that have grown more 
massive than M ~ 10 13 Mq and use an JV-body simulation 
to follow their accretion onto larger dark matter halos with 
z — mass comparable to observed rich groups and clusters 
(M > 10 14 Mq). Comparison with the observed multiplicity 
function of LRGs at z ~ 0.3 implies that accreted LRGs 
must merge on timescales comparable to those predicted by 
Chandrasekhar's formula (~ 2 Gyr). While this may not be 
surprising, the relative flood of massive halos onto more mas- 
sive halos implies that a substantial number of LRGs have 
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Figure 1. Fraction of halos at z = with M > 10 14 Mq, /, that 
have absorbed -/V sa t halos with mass > 10 13 Mq . The thick solid 
line is for halos accreted at all epochs, while the dashed (dotted) 
line indicates only those halos accreted after (before) z = 0.8, 
i.e. when the Universe was half its present age for our assumed 
cosmology. No distinction is made here between halos that have 
dissolved and halos that remain as bound satellites. On average, 
3.2 halos with mass > 10 13 Mq have accreted onto these more 
massive halos by z = 0. 



disrupted over the history of the Universe. Though these 
numbers may at first glance appear large (on average 2 — 3 
disrupted LRGs per z = halo with M > 10 14 M Q ) we 
show that the total stellar mass brought in by these accreted 
LRGs is consistent with the observed stellar mass in clus- 
ters so long as one counts both observed massive galaxies 
and the observed intracluster light. 

The following sections describe in more detail the salient 
accretion properties of massive dark matter halos (SJ2), the 
inferred merging timescale of LRGs, if LRGs correspond 
to massive halos ($3}, and the implied total stellar mass 
brought into massive z ~ dark matter halos by these 
accreted LRGs We conclude in iJS] Throughout we 

assume a fiat ACDM cosmology with (f2 m , fijy, h, (Ts) = 
(0.25, 0.75, 0.72, 0.8), and use a virial mass definition, M200, 
corresponding to the mass contained within a region that 
has mean density e qual to 200 x the critical density (see e.g. 
lEvrard et ai1l2007l ). 



2 THE ACCRETION HISTORY OF MASSIVE 
DARK MATTER HALOS 

A robust expectation of a Universe dominated by cold dark 
matter is the hierarchical growth of structure, and in partic- 
ular the growth of dark matter halos via the accumulation of 
smaller halos. An illustrative example of the accretion his- 
tory of dark matter halos is shown in Figure Q] There we 
plot the multiplicity function for halos more massive than 
10 Mm, i.e. the distribution of the number of halos with 



mass greater than 10 14 Mq (see the appendix for details re- 
garding the simulation used to compile this information). 
We refer to these more massive halos as "hosts" through- 
out. There are 2339 such hosts in our simulation at z — 0, 
corresponding to a number density of ~ 2 x 10~ 5 Mpc -3 . 
Note that these distributions are not symmetric. On average, 
halos more massive than 10 14 Mq have been bombarded by 
3.2 halos with mass > 10 13 Mq over a Hubble time. Here we 
do not distinguish between halos that were accreted directly 
onto the host halo and those that were accreted onto an in- 
termediate halo that later accreted onto the host, although 
such a distinction will be utilized in the following sections. 
The accretion of such massive halos is roughly equally im- 
portant both at low and high redshift: on average two such 
halos have been accreted at z < 0.8 (the Universe was about 
half its present age at 2 = 0.8 for the cosmology assumed 
herein) . 

Halos more massive than ~ 10 13 Mq are expected 
to contain at least on e massive galaxy at their cen- 
ter (|Zehavi et al ] l2005lj ) even at mod erate redshifts (e.g. 
lYan et al.ll2003l , |2004| ; ICoil et al.|[200f3 ). From Figure [T] we 
are lead to the conclusion that, in the absence of merg- 
ers, observed clusters with M > 10 14 Mq should contain 
on average 3.2 massive galaxies (and certainly more if ac- 
creted halos of lower masses also contain massive galaxies), 
with a significant tail toward much larger numbers. How- 
ever, reproduc ing the observed clu stering of massive galax- 
ies at z ~ |Zehavi et al.l [20 05a) would require closer to 
1.2 galaxies in s uch halos in our simulation, in agreemen t 
with other work l|Masiedi et al.ll200rj ; Ikulkarni et al.ll2007r i. 
While these statements are only qualitative, they will be 
confirmed in the more quantitative discussion that follows. 
In order to reconcile the accretion properties of halos with 
observations, we are thus lead to consider the fate of these 
massive halos and the galaxies within them. 



3 THE MERGING TIMESCALE OF LRGS 

LRGs are massive galaxies with very little ongoing star- 
formation; they thus constitute the tip of the red se- 
quence. They have uniform spectral energy distributions 
marked by numerous features and hence their redshifts 
are relatively straightforwa rd to estimate photometrically 
|Padmanabhan et al.l [20051 . redshift uncertainties are Sz ~ 
0.03). Modeling of their spectral energy distributions has 
lead to the conclusion th at these galaxies formed the bulk 
of their stars at z > 2 (e.g . iTrager et alj|200ol ; Ijimenez et al.l 
2006; Thomas et al.|[2005l ), and hence are expected to evolve 
largely dissipationlessly at z < 1. Their clustering strength 
is large, suggesting that they live in m assive dark matter 



halos M > 10 13 Ma, JZehavi et al.ll2005al ) 



M > 10 13 Mq that have accreted onto halos with z 







Recently, iHo et all (|2007h has measured the multiplic- 
ity function of LRGs extracted from the Sloan Digital Sky 
Survey (SDSS; lAdelman-McCarthv et ai]|2006ft for 43 clus- 
ters over the redshift range 0.2 < z < oTsET 

Cluster virial 

masses were derived from ROSAT X-ray data and range 



The primary spectral feature used to measure photometric red- 
shifts of LRGs is the 4000A break; at z < 0.2 this feature moves 
out of the SDSS bandpass filters. Hence our sample is restricted 
to z > 0.2. Although Ho et al.'s sample extends to 2 ~ 0.6, for 
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from 10 14 ' 1 < M200/M© < 10 14 9 . The average stellar mass 
of the LRGs in this sample is M sta r = 10 n ' 6 MQ, as de- 
termined from a color-based st ellar mass estimator (for a 
Chabrier IMF; iBell et al.ll200ot ). These clu sters contain on 
average 2.5 LRGs. The reader is referred to lHo et all (|2007T ) 
for further details regarding these observations. 

In this sample there are approximately five clusters that 
contain no LRGifl When plotting the observed multiplicity 
function below we include both the one reported in Ho et al. 
and one where clusters with JV < 1 are artificially assigned 
N = 1. This is done to afford a more robust comparison to 
our simple model (see below) where we assume that each 
cluster halo contains at least one LRG at the center. As 
discussed below, our conclusions are insenstive to this dis- 
tinction. 

The dark matter halo accretion history of massive halos 
(e.g. Figure [T]) is closely related to the LRG multiplicity 
function. The former can be converted into the latter if one 
knows both the minimum halo mass (measured at the epoch 
of accretion) associated with accreted LRGs, M m j n , and the 
average time it takes for LRGs to merge and/or disrupt 
once accreted. Below we argue for reasonable values of M m i n 
and then attempt to directly constrain the average LRG 
merging timescale. We parameterize the probability that an 
LRG will have merged by a time £ acc since accretion onto 
the host halo via: 



P — 1 

J merge — - 1 - 



(1) 



where r is the merging timescale. The number of LRGs pre- 
dicted by this simple model is then 



AW = 1 + ^2 e 



1/1 



(2) 



where the first term counts one LRG at the center of the 
host halo and the second term counts those satellites with 
accretion epoch mass > M m j n that have not merged. For the 
purposes of generating a multiplicity function we round the 
second term to the nearest integer. Note that in generating 
a multiplicity function we do not have to make a distinction 
between accretion events that did or did not occur within the 
main host halo. This distinction will only become relevant 
when discussing the merger rates of LRGs. 

The minimum LRG stellar mass in the observed sam- 
ple is 1O 113 M . A minimum halo mass associated with 
LRGs can be estimated by assuming the universal baryon 
fraction ft = 0.17 and an efficiency factor r\ of convert- 
ing baryons into stars. While this factor is only poorly con- 
strained even at low redshift, values of order ~ 0.1 are likely 
reasonable for these mas sive galaxies l|Hoekstra et al 1 l2005l ; 
iMandelbaum et al.ll2006h and imply a minimum LRG host- 
ing halo mass of M m i n ~ 10 13 M@ . Note that the minimum 
halo mass for hosting LRGs is here used to identify those 



our purposes we truncate it at z = 0.5 to limit the amount of 
possible evolution within the sample. 

2 We say approximately because Ho et al. statistically remove 
interlopers based on photometric redshift uncertainties and hence 
clusters contain a non-integer number of LRGs. 

3 Throughout we use the words "merge" and "disrupt" inter- 
changeably since our analysis does not distinguish between these 
two possibilities. 



T = 2.0 Gyr 




x = 3.0 Gyr 




T = 5.0 Gyr 




Figure 2. The multiplicity function of observed LRGs [dashed 
line joining points) and accreted dark matter halos for two ac- 
creted mass thresholds: M > 10 13 Mq (solid line) and M > 
2 X 1O 13 M0 (dotted line). The observed curves diverge at JV < 1 
due to whether or not we artifically assign an LRG at the cen- 
ter of observed clusters with N < 1 (see text for details). Top 
Four Panels: Accreted halos are assumed to have disrupted af- 
ter a time t, shown in the upper right corner of each panel. It 
is apparent that if LRGs can be identified with halos of mass 
> 10 13 Mq then they must on average merge within ~ 2 Gyr. 
Bottom Panel: Accreted halos are assumed to have merged after 
a dynamical friction timescale (see Equation [3} . This panel also 
includes halos more massive than M > 5 X 1O 12 M0 (dot-dashed 
line) for comparison. 



accreted halos that are likely to host an LRG, with a signif- 
icant fraction of the accretion occurring at 2 ~ 1 or higher. 
Thus this minimum mass will likely not directly correspond 
to the minimum halo mass hosting LRGs at 2 ~ 0, since 
the LRGs that have survived to the present epoch will have 
accreted much more dark matter, resulting in a larger value 
for M m i n at the present epoch. In our simulations halos to- 
day are on average five times more massive than they were 
when they first crossed M m i n , making our est imate consis- 
tent with that inferred from 2 ~ clustering |Zehavi et al.l 
l2005al ; lKuikarni et al.ll2007h . As we discuss in i|4] M m i n much 
larger than 10 13 Mq would require unreasonably long dy- 
namical friction times and M m i n either much larger or much 
smaller would be in conflict with stellar mass estimates in 
clusters. 

Figure [2] plots the resulting LRG multiplicity function 
both for LRGs in accreted halos more massive than 10 13 M 
(solid lines) and observations (dashed lines). We also include 
predictions for LRGs associated with halos twice as mas- 
sive as our fiducial minimum mass (dotted lines) in order 
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$ 1.0 




Figure 3. Number of accreted halos with M > 10 13 Mq that 
have merged with the host halo by redshift 2, per z = host 
halo with M > 10 14 Mq. The merging time is computed in two 
different ways: a constant timescale of either or 2 Gyr (dashed 
and solid lines) and a timescale determined by the Chandrasckhar 
dynamical friction formula (dotted line). 



to illustrate the sensitivity to our assumed LRG halo mass 
threshold. Each panel in Figure [5] is the multiplicity func- 
tion for a different merger timescale. The top four panels 
assume that the merger timescale, r, is constant. It is ap- 
parent from these panels that if M m i n ~ 10 13 Mq then LRGs 
must merge on a characteristic timescale of ~ 2 Gyr. This 
timescale implies an average number of LRGs per cluster of 
2.5, satisfyingly close to the observed value of 2.6. 

Note that in order to compare to the observations we 
have weighted host halos at z = 0.3 (of which there are 460 
in our simulation within the observed mass range) in such 
a way as to reproduce the mass distribution of the observed 
clusters. It is the combination of these two effects (higher 
redshift and different mass distribution) that does not allow 
a direct comparison between Figures [1] and [2] 

If LRGs never merged, there would be on average 5.8 
LRGs per z — 0.3 host halo (averaged over the observed 
distribution of halo masses) . Comparing this number to the 
observed 2.5 LRGs per cluster highlights the importance and 
prevalence of LRG mergers. 

The bottom panel in Figure [5] assumes that the merger 
timescale is equal to th e dynamical friction timescale 
|Binnev fc Tremaind[l987l ): 



toF = 0.1 tjr 



Mh/M s 



\n(l + M h /M a 



(3) 



where Mh and M s are the host and satellite masses, z is 
the redshift, tn is the Hubble time and all quantities are 
measured at the epoch of accretion. The average mass ratio 
at accretion is ~ 6 for our sample. The pre-factor, 0.1 tn, is 
the characteristic time for a halo with mean density O(10 2 ) 
times the critical density. It is important to note that £df 
gets shorter both at higher redshift and for merger mass 
ratios closer to unity. In this lower panel we have additionally 
included results for halos of mass M > 5 x 10 12 M Q (dot- 
dashed line) for comparison. 

This comparison with simple dynamical friction esti- 



mates provides a satisfying cross check to the results in the 
upper panels. In particular, for our fiducial minimum LRG 
halo mass of 10 13 Mq , the dynamical friction timescale av- 
eraged over all the accreted halos is 2.4 Gyr (median time 
is 1.8 Gyr — the distribution is highly asymmetric), which 
is quite similar to the constant merger timescale that best 
matches the observed multiplicity function (~ 2 Gyr). From 
Figure [2] it is apparent that the simple DF timescale would 
not have reproduced the observed LRG multiplicity function 
if the minimum halos mass capable of hosting LRGs were 
substantially more or less massive than 10 13 Mq. The im- 
plication here is clear: if M m j n is in fact considerably larger 
or smaller than 10 13 Mq then simple DF arguments do not 
apply to the LRG population. 

In fact, it is not at all clear that Equation [3] should ap- 
ply here or in general to the dynamical evolution of satellite 
galaxies, as it is strictly valid for a point mass moving in 
an infinite, uniform background density field. Indeed, much 
work has gone into bot h testing the validity of Equation 
wit h simulations (e.g. [White! 19831: Ivan den Bosch et al.l 



19991 : IVelazouez fc Whitl Il999l; iRead et all |200d) and de- 
veloping extensio ns to it (e.g. [Tremaine fc Weinberg! Il984l : 
IColpi et al.l 1 1999T ) . including numerically following the evo- 
lution of the sat ellite orbit, includin g the mass lo s s due 



to tidal forces dBenson et all l2002l: |Taffoni et all 120031 : 
iTavlor fc Babull |2004 IZentner et'aiTl2005l '). This body of 
work has shown that Equation [3] is at best a crude approx- 
imation to the realistic, time-dependent problem. For these 
and other reasons it is quite surprising, if not entirely co- 
incidental, that the classical Chandrasekhar DF timescale 
adequately captures the merging timescale of LRGs. 

As can be seen from the upper four panels in Fig- 
ure [2] there is a degeneracy between the merger timescale 
and the minimum halo mass associated with LRGs in the 
sense that a larger M m i n coupled to a larger timescale can 
produce roughly the same multiplicity function. Thus, if 
one thought that LRGs lived in more massive halos than 
what we have assumed here, then one would infer a longer 
merging timescale for LRGs. However, this is exactly oppo- 
site to what one would infer from dynamical friction argu- 
ments since tnF oc M s _1 . Furthermore, increasing M m i n to 
5 x 10 13 Mq would result in far too few LRGs in massive 
halos compared to observations, even if r = 00; in this case 
the average number of LRGs per halo would be 1.6. As we 
describe in ^4j M m i n is further constrained by observations 
of the stellar light in massive halos. 

These merger timescales can easily be cast into a dis- 
cussion of LRG merger rates. For this discussion we consider 
the full host halo population at z — 0, rather than the pop- 
ulation at z = 0.3 meant to coincide with the data from 
Ho et al., in order to draw more general conclusions about 
LRG mergers. Figure [3] plots the cumulative distribution of 
LRGs as a function of redshift, per host halo. The 



figure includes constant merger timescales of and 2 Gyr 
and a timescale set by dynamical friction. The Gyr case 
can equivalently be thought of as the distribution of ac- 
creted LRGs, since in this case the accretion and merging 
epochs are coincident. In this figure we only count LRGs 
that merge within the main progenitor of the z — host 
halo. This figure is thus not directly comparable to Figure 
[1] In other words, if a halo merges within a halo that itself 
later mergers with the host halo then it is not counted here. 



© 0000 RAS, MNRAS 000, 000-000 



Constraints on the Merging Timescale of Luminous Red Galaxies, Or, Where Do All the Halos Go? 




M host [ivy 

Figure 4. The total stellar-to-virial mass ratio as a function 
of z = halo virial mass. Only galaxies more massive than 
Mstar = 1O 113 M0 and the ICL are included in the stellar mass 
budget. Observations from iLin et ah! J2004I . solid circles) have 
been converted to total stellar masses, and are compared to to- 
tal stellar masses estimated by assuming that every accreted halo 
with mass > 10 13 Mq hosts an LRG with Afgtat = 10 11 ' 6 M Q 
{diamonds). 

This plot is thus meant to capture the number of mergers 
actually occurring within the main progenitor of the host 
halo. 

The dynamical friction timescale is shorter than 2 Gyr 
at high redshift and longer than 2 Gyr at low redshift; this 
results in a more gradual increase in the merger rate per unit 
redshift compared to a constant merger time of 2 Gyr. Since 
the constant and dynamical friction timescales are different 
at redshifts both greater and less than 0.3, comparisons to 
the multiplicity function at different epochs can in princi- 
ple rule out either (or both) of these timescales. Both the 
constant 2 Gyr timescale and that determined by dynamical 
friction imply that 2 — 2.5 LRGs have merged with the host 
halo by z — 0. Moreover, the figure indicates that a sub- 
stantial number of LRGs are merging/disrupting at z < 1. 
In the next section we set this in the context of recent ob- 
servational results of the stellar mass budget in groups and 
clusters. 



4 THE TOTAL CLUSTER STELLAR MASS 

The merger rate of LRGs found in the previous section in- 
dicates that a significant number of LRGs must have dis- 
rupted over the history of massive halos. In particular, a 
merger timescale of r = 2 Gyr implies that on average 2 — 3 
LRGs have disrupted within each host halo more massive 
than 10 14 Mq by z = (and approximately 2.1 since z = 1). 
In this section we discuss how reasonable such a disruption 
rate is in light of the total stellar mass observed in clusters 
at z ~ 0. 

The amount of stellar mass brought into massive halos 
by LRGs can be estimated in the following way. We consider 
halos that have a mass at z = greater than 10 14 Mq and 
assign a stellar mass of 10 116 M Q , the mean stellar mass of 



the observed LRG sample, to the center of each. Then, each 
halo that was accreted onto these z — halos is assigned 
the same amount of stellar mass if the halo mass at accre- 
tion is > 10 13 Mq . This exercise is thus meant to count the 
total amount of stellar mass that was at some point associ- 
ated with LRGs. We make no distinction between disrupted 
and non-disrupted LRGs except in one case: if, according 
to our best-guess LRG merger timescale (2 Gyr), an LRG 
halo would have disrupted not in the main progenitor of 
the z — halo but rather in some smaller halo that would 
itself later accrete onto the main progenitor but does not 
merge, then we do not count this LRG in the final stellar 
mass budget. In this case the disrupted LRG contributes to 
the satellite's ICL (i.e. ICL that surrounds the satellite and 
is distinct from the central ICL). As it turns out, only 20% 
of all accreted halos more massive than 10 13 Mq fall into 
this category and including these halos in the stellar mass 
budget does not appreciably change our conclusions. 

We compare to data presented in ILin et al l (|2004h who 
have compiled information on 93 clusters at z < 0.1, includ- 
ing X-ray observations used to derive cluster virial masses 
and luminosities of cluster members derived from 2MASS 
photometry. From these data Lin et al. have estimated the 
luminosity function (LF) of each cluster, assuming that the 
faint-end slope is fixed at a — —1.1. Using their LFs we 
are able to estimate the total luminosity in galaxies brighter 
than Lk = 2.8 x 10 11 Lq which corresponds to the minimum 
luminosity of the Ho et al. LRG sample. This total lumi- 
nosity is converted into stellar mass by assuming a mass- 
to-light ratio of M^^/Lk = 0.72 w hich is appropria te for 
red galaxies wi th a Chabrier IM F (|Bell et al) boost ) . The 
LFs reported in lLin et al.l l|2004l ) do not include the bright- 
est cluster galaxy (BCG); we thus add these in separately. 
Finally, we have assumed that each cluster contains intra- 
cluster light (ICL) with a stellar mass equal to the mass 
of the BCG identified by Lin et al. (i.e. Licl = £bcg)- 
This amount of light ass ociated with the ICL is consistent 
with recent observations (|Zibetti et al . 2005; Gonzalez et all 
2005). In this "total" cluster luminosity we do not include 
any possible ICL associated with satellites; it is for this rea- 
son that we did not include LRGs that disrupted in halos 
which themselves later merged with the host. 

Figur e U presents a comparison between the data from 
I Lin et all l|2004h and the stellar mass associated with ac- 
creted halos. The agreement is encouraging. Note that vary- 
ing any one of our assumptions can change the results from 
both the data and our model; the important point to take 
away from this comparison is that the influx of massive 
galaxies embedded within accreted halos appears to roughly 
agree with the total stellar mass within observed clusters 
at z ~ 0. This provides further support to our identifica- 
tion of halos with mass > 10 13 Mq as being host to LRGs 
and suggests that disrupted LRGs deposit their stars into a 
combination of the central galaxy and ICL. Increasing or de- 
creasing M m i n by a factor of two would result in substantial 
disagreement with the observations shown in Figure [4] This 
is due to the fact that the number of accreted halos does 
not scale linearly with the accreted halo mass, and provides 
further support for our choice of M m t n = 10 13 Mq. 

The simple model presented here also provides a 
straightforward means for understanding the observed trend 
of decreased scatter in M star /Mhost with increasing Mhost. 
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This arises because the number of accreted halos with mass 
is a weak function of Mhost ■ This is in contrast to 
the observed num ber of satellites, which appears to be closer 
to linear in M hoa t (|Lin et alJl2004|Popesso et alj|2007l ). The 
difference implies that fractionally more halos/LRGs are 
merging in lower mass halos compared to higher mass halos. 

Interestingly, we find that it is not uncommon for z = 
host halos to contain disrupted LRGs that did not disrupt 
within the host halo (and hence were not counted in the 
above figure) but rather disrupted in a smaller halo that 
later accreted onto the host (and remained as a satellite to 
z = 0). If these disrupted LRGs are depositing some fraction 
of their stars into ICL, then this suggests that there could 
be a significant amount of ICL that is not centered on the 
central galaxy but is instead centered on cluster satellites. 
Such a scenario is corroborated by recent observations of 
the Virgo cluster that show significant amount s of ICL sur- 
roun ding several of the most massive satellites (|Mihos et al.l 
l2005h . 



5 DISCUSSION 

The results of the previous sections suggest the following 
picture. If LRGs are associated with halos more massive 
than 10 13 Mq at the time when they are accreted onto more 
massive host halos, then the observed multiplicity function 
of LRGs at 2 ~ 0.3 implies that LRGs must merge and/or 
disrupt on timescales of ~ 2 Gyr. Such a merger rate im- 
plies that 2 — 3 such LRGs have disrupted in halos more 
massive than 10 14 Mq by z = 0. This merger timescale is 
consistent with classical dynamical friction arguments and 
suggests that a rather simplistic dynamical prescription for 
the evolution of LRGs is applicable when considering en- 
semble averages. 

Moreover, the amount of total stellar mass in clusters 
that was at one point associated with these infalling LRGs 
(ignoring for the moment whether or not this stellar mass 
is locked up in satellite galaxies) is consistent with observa- 
tions when the observed amount of stars in the intracluster 
light (ICL) is accounted for. This in turn suggests that the 
disrupting LRGs are depositing their stars into a combi- 
nation of the ICL and central galaxy, which is consistent 
with previous modeling ( Monaco ct al. 2006; Mu rante et alj 
120071 ; iPurcell et al.ll20OT ; IConrov et al.ll2007l ). Finally, there 
appears to be a significant number of LRGs that have dis- 
rupted within halos that only later accreted onto (but did 
not merge with) what would become the z — host halo. 
This suggests that there could be a significant amount of 
ICL surrounding cluster satellites, in addition to what is 
known to be associated with the central galaxy. 

It has been historically challenging to constrain the 
merger rate of galaxies. Previous stu dies have relied on ei- 
ther morphologica l disturbances (e .g. IConselice et al.ll2003l ; 
van Dokku m 2005; Bell et ail l2006) or close pair counts (e.g. 
Masiedi et al.l 12006) as probes of the merger rate of mas- 



and z = 0.2, the model presented herein implies an LRG- 
LRG merger rate of (1.0 - 1.3) x 10 4 Gyr -1 Gpc -3 , de- 
pending on whether the constant 2 Gyr or dynamical fric- 
tion t imescale is used. The agreement with Masiedi et al.l 
(2006) is encouraging, especially given the (different) uncer- 
tainties in both approaches. These rates are also consistent 
with curren t predictions from cosmological hydrodynamic 
simulations l|Maller et alj|2006h . 

Many studies have attempted to constrain the stellar 
mass growth of massive gal axies from their infe rred merger 
rates. However, as argued in IConrov et ail (|2007h and herein, 
the merging of massive galaxies will often not correspond to 
significant growth of the resulting galaxy because a substan- 
tial amount of stars can be transfered to the ICL. We hence 
caution against using merger rates to constrain the stellar 
mass growth of galaxiesQ In fact, significant growth of the 
ICL via merging at late times provides a means for reconcil- 
ing two apparently contradictory facts: one the one hand, ob- 
servations at z < 1 indica te that central mass i ve red galaxies 
grow little in mass (e.g. iBrown et al.l 120071 ; iFontana et al.l 
120061 ; iBundv et all 120061 ; IWake et al.l l2006h . while on the 
other hand, merging/disruption of galaxies within groups 
and clusters at la te times appears relatively common (e.g. 
I White et afl 120071 ). 



White et all (120071 ) outlined an approach for measuring 



sive galaxies. Unfortunately, both methods are rather in- 
direct since the connection between either morphological 
disturbances or close pair counts and merger rates is un- 
certain. _rhe_jnost_recent_Jxiferred LRG-LRG merger rate 
is from iMasiedi et ail l|2006h who find a rate of 0.6 x 10 4 
Gyr -1 Gpc -3 . Averaging over all halos between z — 0.5 



the merging rate of massive galaxies similar to the one pre- 
sented herein. Using the observed evolution in the cluster- 
ing of massive galaxies, these authors concluded that ~ 1/3 
of massive satellites merge/disrupt between z ~ 0.9 and 
z ~ 0.5. In the present work we find roughly 50% of massive 
satellites have disrupted over similar epochs. Our fraction 
is slightly highe r because we have f ocused on more massive 
galaxies than in IWhite et al.l(|2007f ). The more general con- 
clusion from these two studies is, however, robust — the pop- 
ulation of massive galaxies experiences significant amounts 
of merging/disruption, even at z < 1. 

There is an important implication of considering the 
evolution of galaxies within the context of the hierarchical 
growth of halos. At first glance, the lack of evolution in the 
observed correlation function of massive galaxies and their 
halo occupation distribution at z < 1 suggests that massive 
galaxies do not disrupt or merge over this epoch. However, 
these galaxies are embedded within dark matter halos that 
are continually merging and accreting new galaxies, which 
instead suggests that massive galaxies must merge in or- 
der that these observed quantities not evolve appreciably at 
late times. This statement is further corroborated by dissi- 
pationless simulations which show explicitly that the aver- 
age number of subhalos within host halos does not evolve 
appreciably at z < 1 because the accreti on and disruption 
rate of subhalos ar e approximately equal (|Reed et al. I l2005l ; 
IConrov et al.|[20 06). If satellite galaxies reside within these 
subhalos then the observed non-evolution of the clustering 
and halo occupation of massive galaxies at 2 < 1 is in fact 
consistent with significant amounts of merging at late times. 

Our results highlight the power of using purely dissi- 
pationless simulations coupled to simple relations between 



4 This issue is intimately related to the way in which one counts 
galaxy light. Of course the combined light of both the central 
galaxy and its ICL will increase after a merger event. 
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galaxies and dark matter to infer the evolution of galax- 
ies and their relation to the underlying dark matter with 
time. The approach outlined herein can easily be extended 
to other datasets to provide additional constraints on the 
merger rate of galaxies. 
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APPENDIX A: THE SIMULATION, HALO 
CATALOG, AND MERGER TREES 

We use a high resolution simulation of a ACDM cosmology 
(Q. M = 0.25 = 1 - O a , «s = 0.043, h = 0.72, n s = 0.97 and 
erg = 0.8). The linear theory power spectrum is computed 
by evolution of the coupled Einstei n, fluid and Boltzman n 
equations using the code described in l White fc Scottl (1 19961). 



This code agr ees well with CMBfast ( Seliak fc Zaldarriagal 
Il996h . see e.g. ISeliak et alj (|2003f ). The simulation employs 
1024 3 particles of mass 8 x 10 9 /i _1 M in a periodic c ube of 
side 500 7i _1 Mpc using a TreePM code l | White! 120021 1. The 
Plummer equivalent softening is 18/i _1 kpc (comoving). 

The phase space data for the particles exists at 50 out- 
puts, spaced equally in conformal time between z ~ 3 and 
2 = 0. For each output we generate a cat alog of halos using 
the Friends-of- Friends (FoF) algorithm l)Davis et al.lll985l l 
with a linking length of 0.168x the mean inter-particle spac- 
ing. This procedure partitions the particles into equivalence 
classes, by linking together all particle pairs separated by 
less than a distance b. The halos correspond roughly to parti- 
cles with p > 3/(27r6 3 ) ~ 100 times the background density. 
For each halo we compute a number of properties, including 
the mass M200 interior to r2oo within which the mean den- 
sity is 200 x the critic al density. M200 is c omputed from a 
fit of an NFW profile (jNavarro et al.lll997i l to the particles 
in the FoF group. 

Merger trees are computed from the set of halo catalogs 
by identifying for each halo a "child" at a later time. The 
child is defined as that halo which contains, at the later time 
step, more than half of the particles in the parent halo at 
the earlier time step (weighting each particle equally). For 
the purposes of tracking halos this simple linkage between 
outputs suffices (note that we do not attempt to track sub- 
halos within larger halos, which generally requires greater 
sophistication). From the merger trees it is straightforward 
to compute the time when a halo 'falls in' to a larger halo, 
the number and masses of the progenitors etc. 



